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PROGRAM DK8CRIPT1QN GUIDE 


#• 

F 

A. IDENTIFICATION 

u rogram Name - Erythropoiesls Regulatory Model 

Programmer and 

Bioengineer Contact - D. J. Grounds or J. I. Leonard, GE/T8SD, 

P. O. Box 58408, Houston, Texas, 77058, 


Date of Issue - March 10, 1978 

B. GENERAL DESCRIPTION 

A mat he iatlcal model and digital computer simulation of the erythropoietic 
control system were developed to simulate long term dynamic changes in red 
cell mass as a result of various environmental and physiological stimuli, 
including those associated with bed rest and weightless spaceflight. The 
purpose of the model is to provide a method to analyze some of the events that 
could account for the decrease in red cell mass observed in crewmen return- 
ing from space missions. 

C. USAGE AND RESTRICTIONS 
Machine, Operating System, 

and Compiler Required - Univac 1110, EXEC 8, FORTRAN V 

Peripheral Equipment Required - ' Time-Sharing Terminal, Tektronix 4010 

Graphics Terminal 

D. PARTICU 1AR DESCRIPTION 

The model is based on the premise that erythrocyte production is governed 
by the balance between oxygen supply and demand at a renal sensing site 
(Figure 1). Oxygen supply is taken to be a function of arterial oxygen tension, 
mean corpuscular hemoglobin concentration, oxy-hemoglobin carrying capacity, 
hematocrit, and blood flow. Erythrocyte destruction is based on the law of 
mass action. The instantaneous hematocrit value is derived by integrating 
changes in production and destruction rates and accounting for the degree of 
plasma dilution. 
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* IMPLEMENTED IN VERSION I MODEL 
O IMPLEMENTED IN VERSION II (CURRENT) MODEL 


FIGURE 1 

THE IMPROVED ERYTHROPOIESIS MODE I 
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The* moat significant improvement in the model is the explicit formulation 
of separate elements representing erythropoietin production and red cell pro- 
duction. Other miNlificat ions include hone marrow timc-dclavs, capability to 
shift oxyhemoglobin affinity and an algorithm for entering experimental data 
as time-varying driving functions. These changes result in a significantly 
Improved and more flexible model as demonstrated by consistently more 
realistic simulations. 

A more detailed description of this model and the results of a validation study 
may be found in References 1-3. 

The various types of simulations which can be performed with this model are 
illustrated in Table 1 

In addition to these input parameters which can be changed to simulate experi- 
mental or clinical conditions, the model can be driven by experimental data 
for plasma volume or hematocrit. Data from experimental studies are avail- 
able to the user upon demand. These studies will be listed under description 
of model input and the references appear in the Reference section. 


E. DESCRIPTION OF INPUT 

(1) Access procedures to the Univac 1110 have been documented in TIH 741- 
MED-3059 "Simplified user instructions for physiological models using 
the Univac 1110 in Demand Mode." 

An example of input for the model is given in Appendix 1. This example 
contains all necessary responses after obtaining a carrier signal on the 
terminal. 

(2) Execution of the program begins by typing (uADD THERM. LOADNEW. 

This instruction loads the appropriate programs into the user's file, 
compiles them, and collects them for execution. This process involves 
several lines of print before the program logins execution. 

Initial conditions for the program are established by reading in a set of base- 
line values for the model parameters. The parameter names, definitions, 
and initial values and units are given in Appendix II. The mnemonics given in this 
glossary are the ones which must be used for specifying changes in input values 
or the output list, (see below). 
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(3) The user in given an option to use experimental datn as Input. If this 
option is exercised, the following choices are offered. 

1 Bed Best Plasma Volume (Morse, 1969) 

2 Bed Rest Plasma Volume (I<li'ali/ed) 

3 Bed Best Plasma Volume (BCM -Johnson, 1970) 

4 Bed Best HCT (BCM- Johnson, 1976) 

5 Altitude Hypoxia Plasma Volume (Buderer, 197.1) 

6 Skylah 3 Plasma Volume 

7 Skylah 3 HCT (Peripheral) 

8 Skylah 4 Plasma Volume 

9 Skylah 4 HCT (Peripheral) 

10 Input Plasma Volume 

11 Input Peripheral Hematocrit 

Upon selecting a study, the stored experimental data are written out. 

(4) The user is then asked to specify the number of days for the run and 
the print Interval. 

The question INPUT WANTED CHANGED gives the user the opportunity 
to change any of the 07 model parameters listed in Appendix II. Some of 
the most meaningful parameters to change upon input are given in Table 1. 

TABLE 1 

INPUT PARAMETERS 
Pefinit Ion 


time limit of run (days) 
print interval (days) 

Integration step size (min) 

Initial red eell volume (liters) 
plasma volume (liters) 
blood flow (I iters/m in) 

multiplier for changing basal oxygen consumption of 
renal tissue (ratio, to normal) 

mean corpuscular hemoglobin concentration 
(gm Hb/llter blood) 


Mnemonic 

TMAX 

TPNT 

H 

RC VO 
PV 
BE 
AOM 

MCHC 
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TABLE 1 (Continued) 
IN PIT PARAMETERS 


Mnemonic 


Definition 

CHB02 

m 

maximum oxygen carrying capacity of hemoglobin 
(ml 0.,/gm lib) 

RKC 

m 

clearance time constant for red cell destruction 
(min -1 ) 

P02ART 

m 

arterial blood pO i( (mm Hg) 

Cl 

m 

sensitivity or gain of erythropoietin releasing system 
to tissue pO 

m 

G2 

■ 

sensitivity or gain of bone marrow to ilasinn erythro- 
poietin concentration 

EH L 

- 

erythropoietin half life (minutes) 

ENIF 

c 

ervthro|x)letin infusion (x normal volume) 

RCVLO 

s 

volume labeled red cells infused (liters) 


F. DESCRIPTION OF Ol'TPUT 

The question Ot’TPl'T WANTED or SAME gives the user the option to replace 
the output list with any specified parameter in the glossary. The default list is 
given in the example developed in Appendix I, 

The model is set up to give graphieal or tabular output on demand. If tabular 
output is selected, the user is guven the opportunity to plot the data in the 
output list after the tabular output is complete. 

The graphics option permits the user to plot the model results using a 
TEKTRONIX -1010. Six ((») parameters may lie selected from the eight selected 
output parameters They are plotted against time in days and the loc (LOC) 
will break up the page into (2-6) separate plots. 

The user is also required to specify the upper and lower coordinate limits 
(HIGH , LOW). A sample run illustrating this procedure is included in Appen- 
dix I 
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G. INTERNAL CHECKS AND EXITS 

The model will simulate desired experiment until time equals TMAX. The 
user will have an opportunity to: 

a) continue the run by typing MORE 

b) begin a new run by typing NEW 

c) exit from execution mode by typing STOP 

If you elect to continue the run (MORE) the present values of the input para- 
meters will be printed and you may change any of them. The very least 
change that is required is to increase TMAX from its current value. Time 
is reset back to zero only when a NEW run is l>egun. 

A normal stop exit occurs by typing "STOP" In response to "WISH TO CON- 
tlnue?" 


H. INDEPENDENT SUBROUTINES 


RE DC EL 

INPUT 

TISSUE 

BLOOD 

OUTPUT 

INTGRL 

GETSP 

PLOT 

XPER 

DATA 1-DATA 9 
BDAT 


Main Program 

Initialization and conversational input 
Computes tissue oxygen tension 
Computes red cell mass 
Conversational output 
Integration algorithm 
Oxv-hemoglobin saturation algorithm 
Graphics algorithms 
Reads experimental data 
Experimental Data 
Initialization Data 


I. 


SAM PLE RUNS 


See Appendix I . 


J . C OMPUTER PROGRAM LISTING 


* 

a 

i 


See Appendix III. 


7 


REFERENCES 


1. Leonard, J. I. "Study Report on Improvements and Validation of the Erythro- 

potesis Control Model for Red Heat Stimulation," Technical Information 
Release, Till 782-LSP-7012, General Electric Company, Houston, Texas, 

1 1*77. 

2. Leonard, J. I. , "Study Report - The Application of Systems Analysis and 

Mathematical M«>del8 to the Study of Erythropolesis During Space Flight, " 
Technical Information Release, TIR 741-MED-4012, General Electric 
Company, Houston, Texas, 1974. 

3. Leonard, J, I. "Validation of a Model for Investigating Red Cell Mass Changes 

During Weightlessness," Technical Information Release, Till 782- MED- 
6001, General Electric Company, Houston. Texas, 1976. 

4. Morse, B. S. "Erythrokinetlc Changes in Man Associated with Bed Rest. " 

Lectures in Aerospace Medicine, 6th Series, School of Aerospace Medi- 
cine, Brooks AFB, Texas, pp, 240-254, 1967. 

5. Buderer, M. C. and N. Pace. "Hemopoiesis in the Pig-Tailed Monkey Macaca 

Ncmcstrlna During Chronic Altitude Exposure," Amer. J. Phvsiol. 223: 
346-352. 1972. 

6. Johnson, P. C. and T. Driscoll. "Red Cell Mass and Body Volume Changes 

After 28 Days of Bedrest," In: "Report of 28-Day Bedrest Simulation of 
Skylab." Volume II, P.C. Johnson and C. Mitchell (EdB. , NASA Contract 
NAS9-I4578, The Methodist Hospital, Houston, Texas, 1977. 

7. Johnson, P.C. , Driscoll, T. B. , and A. D. LeBlane, "Blood Volume Changes, " 

In: Proceedings Skylab Life Science Symposium, Vol. II, NASA TMX-58154, 
National Aeronautics and Space Administration, Lyndon B. Johnson Space 
Center, Houston, Texas, 77058, pp. 195-505, 1974. 


APPENDIX I 
EXAMPLE RUN 
SIMULATION OF BED REST 


T XL 0 06 

ENTER USEPID'PrtiSwORDJ 

> 


♦DESTROY USERID'PA. SWORD ENTRY 

•UN1VAC 1100 OPERATING IYSTEM VER. 31.244.214C • RSI > ♦ 
>i'RUH. N DG629* 7007-05 09-C* DB8-G03432. 1 0 
DA T E » 030378 TIME! 150722 

>i>ADD THERM. LQADHEW 

FURFUR 27R1 RL72-8 03 '09 '78 151071 34 
1 SYM 
1 SYM 
1 SYM 
1 SYM 
1 SYM 

♦OUTPUT INTERRUPT# 

99X □ 

1 SYM 
>9XQ r 

150 CARDS READ AS INPUT 


♦ ♦MODEL FOR REGULATION i’R ERYTHROROIESIS## 

EXPERIMENTAL I'AIA WANTED? cY'N> 

>Y 

DO VO'J WISH INSTRUCTIONS <Y/N> 

>Y 

DATA CAN FE SELECTED FROM THE FOLLOWING STUDIES 
1= FEDREST PLASMA VOLUME (HORSE. 1989) 

2* FED PEST PLASMh VOLUME -IDEALIZED) 

3* FED REST PLASMA VOLUME (l CM- JOHNSON. 1976) 

4= FEDPE1T HCT (ECM- JOHNSON* 1976) 

5* ALTITUDE HYPO IA PLASMA VOLUME cBUDERER* 19?3> 
t= ZKYLAB 3 PLASMA VOLUME 
7 = STYLAB 3 HC T •PERIPHERAL) 

8- SKYLAB 4 PLASMA VOLUME 
9* SKYLAB 4 HCT •'PERIPHERAL) 

10 = input plasma volume 
11= INPUT PERIPHERAL HEMATOCRIT 

ENTER EXPERIMENTAL DATA CODE WANTED 

>1 

.00 3. 00 1.U0 2.54 10. 00 

38. 00 2.70 50.00 2. 70 80.00 


2.49 27.00 2.4c 35.00 

3.03 .00 .00 


2 . 48 
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APPENDIX I (Continued) 


HAFD COPY WANTED <Y» N> . . . 

>N 

HO YOU WISH TO PLOT SUFFERED DP TP? *Y/N CR> 

>V » 

<GD 

<VtRL 

HARD COPY WAMTED<Y*N> . . . 

>N 

<b 

GRAPHIC OUTPUT • Y. ft. S' .TIME INTERVALS* STARTX. TTOPX. <A£*3F?.. 0> . . 


> / 4. 

PCV 


Y 

SCALE 

• A4 . SIX* F 4 . U* £FG. 0) 

PLOT <Y»N.S> 

LOC 

Mll5H 

LOW 

• • • 

>n 

EV 


Y 

SCALE 

•H4. Aa»P4. u*£F6. 0) 

PLOT <Y.N. : » 

LOC 

HIGH 

LOW 

• • • 

>N 

HOT 


1 

SCALE 

• A4.ftX»F4. U»£F6. 0> 

plot 'y»n»c > 

LOC 

HIGH 

LOW 

• • • 

>Y 

1 . 

4-*. 



F08TIS 


Y 

scale 

• A4 « 3X • F 4 . U » £F G . 0 > 

plot <y»n*s » 

LOC 

HIGH 

LOW 

• • • 

>Y 

£. 

£4. 

1 *. 


PV 


V 

scale 

' A4 . 8X* F 4 . 0»£FG. 0> 

PlOT<Y.N.S> 

LOC 

HIGH 

LOW 

• • • 

>Y 

1 . 

3. I 

£. 1 


D 


Y 

SCALE 

<.A4» 3X» F 4 . U.£F6. 0) 

PLOT < »'. N* $ > 

LOC 

HIGH 

LOW 

• • • 

P 


Y 

SCALE 

<A4.3X*F4. 0» £F6. 0> 

PLOT <.Yt N» S> 
>N 

LOl 

HIGH 

LOW 

• • • 

ERY 


Y 

SCALE 

• A4»8X*F4. 0. £FG. 0> 

PLOT <Y» N* S> 

lOC 

HIGH 

LOW 

• • • 

>NN 


GRAPHIC OUTPUT <Y»N.S>. TIME INTERVALS* STAPTX*STOPX. 'AE.3F3. 0> . . 

> 

N 

See Figure A-l 

WISH TO CONTINUE? 'ENTER "NEW* "MOPE"* OP "STOP". > 

>STOP 


A-I-2 


APPENDIX I (Continue*!) 


MU. D9T; WHMTE£'*PPlriT STEP 

1 u. • 1 . 

IMPUT Mt=»n TED CmHNGED* EMD 
MhME VHlUE 


>TMftX 1 1 . 

t TMR.C » 1 1.000 CMPMGED FPDM 11 .UU 0 

>£MD 

OUTPUT UHmTEI 1 UP SHME* STOP 
MHME 


> S TOP 

IiQ YOU 1.1HMT GPfrPHlC I MS T £R 1 ’ CP THBULRP OUTPUT? TY/hO 


>N 


CRTS 

*c V 

BV 

HCT 

. 00 

. 

5 . 00 

4 0 . 00 

1.00 

a. 

4.54 

44 . 05 

t* 00 

3 . 00 

4.53 

44 . 09 

3. 00 

2 . 00 

4.53 

44 . 13 

4 . 00 

1.99 

4.53 

44 . 13 

5 . 00 

i . 99 

4.51 

44 . 13 

6 . 00 

l . N 

4.50 

44 . 13 

r.oo 

l.H 

4.43 

44 . 11 

8 . 00 

l . • f 

4 . 4 " 

44 . 09 

?. 00 

l . 9 ? 

4.46 

44 . 06 

10 . 00 

l . 4 - 

4.45 

44 . 03 

11.00 

l . *5 

4.44 

4 - . 97 

PE SHIFT 

-OUT 

HMP 

PETUPM- 


P 03 T 1 S 

PV 

D 

P 

ERY 

30 . 0 U 

3 . 00 

33 . 00 

33 . 00 

1.00 

31 . 4 c 

£.54 

33.00 

£ 1.67 

.91 

31.53 

3.53 

31.99 

£ 0.45 

.83 

31.53 

3.53 

31.97 

19 . 1 3 

.80 

31.53 

3.53 

31.93 

13 . 01 

.30 

31.53 

3.53 

31.38 

17.10 

.80 

31.53 

3.51 

£ 1 • 33 

16 . 39 

.79 

31.53 

3.51 

31.76 

15.34 

.79 

31.53 

3.50 

£ 1.69 

15.41 

.00 

31.53 

£.50 

£ 1 . 63 

15 . 09 

.30 

31.50 

3.49 

£ 1.55 

14.34 

.30 

31 . 4 ? 

3.49 

31.43 

14.66 

.80 




J 1 . 
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FIGURE A-l 

EXAMPLE GRAPHICS OUTPUT 




APPENDIX II 



Mnemonic 

Initial \alue 

i 

T 

0.0 

2 

TMAX 

0.0 

3 

TPNT 

1.0 

4 

H 

1.0 

5 

BF 

1.2 

6 

PV 

3. 0 

7 

RKC 

7. 6389E-06 

8 

EGAIN 

0.0 

9 

P02DMD 

- 

10 

RCVO 

2. 0 

11 

VO20 

0. 032 

12 

BMOTIS 

20. 0 

13 

AOM 

1.0 

14 

P02ART 

95. 

13 

QG2TIS 

C 

16 

P02TIS 

C 

17 

STEP 


18 

BV 

5. 0 

19 

RCV 

C 

20 

HCT 

c 

21 

PROD 

- 

22 

RCD 

- 

23 

DRIVE 

- 

24 

S02ART 

c 

25 

S02VEN 

C 

26 

P02VEN 

c 


Definition 

Time (mins. ) 

Maximum time (days) 

Print Step (days) 

Time Step (reins. ) 

Renal blood flow (L/mln) 

Plasma volume (titers) 

Rate Constant RC Destruction Rate 
(l/mln) 

Not used 

Vf If 

Initial red cell volume (l) 

Basal metabolic rak renal tissue(ml O /min) 

Autonomic stim. effect on tissue C), ; 

Utilization (Normalized) 

PO Arterial Blood (mm Hg) 

O Delivery Rate to Renal Tissue 
(ml/mln) 

O Tension Renal Tissue (mm !’?) 

mt 

Not Used 

Blood Volume (liters) 

Red Cell Volume (liters) 

Hematocrit whole-body 
Not used 

Red Cell Destruction (1) 

Not Used 

0,, Saturation in A-teries (%) 

O l} Saturation in Veins (%) 

Partial Pressure of C) <; in the Veins 
(mm Hg) 


C = Computed values 
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APPENDIX 11 (Continued) 


No. 

Mnemonic 

Initial Value 

Definition 

27 

M02T1S 

20.0 

Metabolic Tissue Rate (ml min) 

28 

P 

22.0 

Hod Cell Production (ml/dny) 

29 

VG2T1S 

320. 0 

Volume of O in Tissue (ml) 
2 

30 

MCHC 

3. 75 

Mean Corpuscular Hq cone, (g m/I) 

31 

CH D()2 

1.34 

Carrying capacity of hemoglobin 
(ml O / kiii Ilg) 

32 

RK 

C 

Diffusion Constant 

33 

I) 

22. 0 

Red Cell Destruction (ml/day) 

34 

DATA 

0.0 

Internal switch 

35 

RCVL 

0 . 0 

lied Cell Volume lain? led (liters) 

3(i 

RCVLO 

0 . 0 

Initial Dose labeled Hod Cell Volume 
Change (liter) 

37 

OPT 

0 . 0 

Internal switch (not used) 

38 

XI 

3. 0 

Values of parameters during o|>ii- 
mlzation Var. 1 

39 

X2 

9.5485 E- 

Values v>f parameters during t>|8 l— 
mlzation Var. 2 

40 

Z 

5760. 

Red Cell Production Delay Time 
Const, (min.) 

41 

P50 

26.73 

Partial pressure O i} (<i 50', saturation 
(mm llg) 

42 

DPG 

37.0 

Dummy variable for P.-„ shift due to 
DPC. 

43 

C.l 

3. 0 

Sensitivity or gain of erythropoietin 
releasing system to tissue pO )T 

44 

G2 

2.0 

Sensitivity or gain of bone marrow to 
plasma erythropoietin cone. 

45 

PO 

0 . 0 

Ilasal red cell production rate 
(normalized) 

46 

PI 

1.0 

Control value of RC Prod, rate (normal 
operating point) (normalized) 

47 

Ell L 

720. 0 

Erythropoietin half-life (min.) 


C Computed Values 
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APPENDIX II (Continued) 


No. 

Mnemonic 

Initiul Value 

48 

ERYO 

1.0 

49 

EPO 

20.09 

50 

ENIF 

0 . 0 

51 

KCELL 

1.0 

52 

DELTI 

0.2 

5.1 

DE LTI1 

1.0 E-06 

54 

EP81 

.01 

55 

E PSI1 

10 E-08 

56 

Dll M4 

50 

57 

DK 


58 

P02N 

1 . 

59 

EP 

1 . 

60 

ED 

1 . 

61 

ERY 

1 . 

62 

RC PN 

1 . 

6.1 

RC P 

. 00001.18 

64 

RC PX 

c 

65 

PRODX 

c 

4 ’ 

PRODZ 

- 

67 

RNVAR 

150. 

68 

FLAG 

0 . 

69 

EX RCM 



Definition 

Initial value of E p (normalized) 

Intercept (changes normal o|iera- 
ting point) 

Erythropoietin Infusion (x normal vol.) 

Ratio of peripheral tICT to whole- 
body HCT 

Parameter, optimization constant 

ft M «« 

Not used 
Not used 
Not used 
Not used 

Tissue pO_ (normalized) 

Erythropoietin production (nor- 
malized) 

Erythropoietin Dost ruct ion 
normalized) 

Erythropoietin Concentration 
(normalized) 

Gain of erythropoietin on red 
cell production (normalized) 

Red cell production (1/m in) 

Red cell production with time dclay(I/min) 

Red cell production with time delay 
(ml/day) 

Not used 

No variables input/output 

Internal switch (RCM data) 

Experimental red cell mass 
(liters) 


C Computed Values 
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